Growth mechanisms of GaN on the O-terminated ZnO(OOOl) surfaces 
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We have investigated the stability of the IML-GaN on the O-polarity ZnO(OOOl) interface struc- 
ture using the first-principles calculation. We have found in our calculated results that the most 
stable structure for the IML-GaN on the O-polarity ZnO(OOOl) interface has the N-polarity. How- 
ever, we have found that the results of the adatom dynamics on the O-terminated ZnO(OOOl) surface 
shows the Ga-polarity. We find that the key to change the polarity of GaN crystal grown on the 
O-terminated ZnO(OOOl) surface is the growth temperature. We have suggested that the opti- 
mized initial growth condition for the growth of the Ga-polarity GaN crystal on the O-terminated 
ZnO(OOOl) surface is under the suitable low temperature and the stoichiometric growth condition. 
Experimental observations show that GaN grown on ZnO(OOOl) by PLD at substrate temperatures 
below 300° C has Ga-polarity, which is quite consistent with the theoretical calculations. 

PACS numbers: 81.15.Fg, 81.15.Aa, 81.05.Ea, 81.10.Aj 
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I. INTRODUCTION 



GaN has attracted much attention as a material 
for short-wavelength optical devices such as blue light- 
emitting diodes (LED) and laser diodes (LD) since it has 
a direct band gap of 3.4 eV.|l|l| Although, GaN films 
are usually grown on sapphire substrates, it is well known 
that these GaN films suffer from formation of a high den- 
sity of structural defects due to the large lattice mismatch 
between GaN and sapphire. ZnO has been regarded as 
one of the most promising substrates because ZnO and 
GaN perfectly share the same crystalline symmetries and 
the lattice mismatches between them are as small as 1.9% 
and 0.4% for the a-axis and c-axis, respectively. In fact, 
Kobayashi et al.Q have recently shown that high qual- 
ity GaN grows in the layer-by-layer mode on ZnO(OOOl) 
even at room temperature (RT) by the use of pulsed laser 
deposition (PLD). They have also shown that the sur- 
face of GaN grown on ZnO (0001) has straight steps with 
a height of one unit cell and large atomically flat ter- 
races. These striking features can be attributed to the 
small lattice mismatch between GaN and ZnO. In spite 
of this success in heteroepitaxial growth, the initial stage 
of the film growth of GaN on ZnO (0001) has not been 
understood very well. Namkoong et al.|J| have recently 
reported that MBE growth of GaN on atomically flat 
O-terminated ZnO (0001) substrates at a substrate tem- 
perature of 600° C leads to formation of mixed polarity 
GaN. However, existence of large atomically flat terraces 
that Kobayashi et al. 3] have reported for RT grown GaN 
implies possibility of growth of Ga-polarity GaN. Fur- 
ther study is necessary to clarify the initial stage of film 
growth and to understand mechanisms that determines 



the polarity of GaN grown on ZnO (0001) substrates. 

First-principle calculation has been proved to be a 
powerful tool for investigation of growth mechanisms of 
group III nitrides. For example, the first-principle cal- 
culations have been successfully applied to investigate 
mechanisms for the Ga-terminated GaN(OOOl) homoepi- 
taxial growth, iQ] Polarity change of GaN by the use 
of the Al-terminated sapphire(OOOl) has also been well 
explained by the first-principles calculation. In this 
study, we discuss mechanisms of the polarity determina- 
tion of GaN on the O-terminated Zn(OOOl) surface with 
the theoretical first-principles calculations of the total en- 
ergy for the IML-GaN on the O-terminated ZnO (0001) 
interfacial structure. 



II. COMPUTATIONAL METHOD 

Our first-principles density-functional total-energy 
calculations are based on the density-functional 
theory (DFT),[i using the Vienna ab initio simu- 
lation package (VAS P).llOl| We used the local-density 
aDProximation(LDA) [l]lTl2 | for the exchange correlation 
and projector augmented wave (PAW) [Hill potentials. 
The Zn and Ga 3d electrons were treated as part of 
the valence band. The cutoff energy for the plane-wave 
basis was 400 eV. The twin boundaries are modeled with 
supercells. We employed supercells containing eight 
atomic layers of ZnO and two layers of GaN, a vacuum 
region 12 A thick. The cation-terminated bottom layer 
of the slab was passivated with fractionally charged 
hydrogen atoms. Two layers on this side of the slab 
were kept fixed to simulate the constrains on the few 
outmost layers coming from the underlying semi-infinite 
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TABLE I; Total energy valuses of IML-GaN on the O- 
terminated ZnO(OOOl) interfacial structure. The energy values 
for the stable structure is set at zero. The unit of the energy 
is eV. 



GaAs/ZnO(0001) structure 


Relative energy (eV/lxl) 


Ga-polarity Wurtzite 


0.93 


Ga-polarity Zincblende 


0.90 


N-polarity Wurtzite 


0.00 


N-polarity Zincblende 


0.08 



bulk, whereas all other atoms were allowed to relax until 
the Hellmann-Fcynman forces became smaller than 0.05 
eV/A. We employed sets of special k points equivalent 
to 4 kpoints within the surface Brillouin zone (BZ). 



III. RESULTS AND DISCUSSION 

We calculated the total energy for the IML-GaN on 
the 0-terminated ZnO(OOOi) surface. We calculated 
four structural models that the polarity has the Ga- 
polarity and N-polarity, and the interfacial structure has 
the Zincblend type and Wurtzite type structure each 
other shown in Fig. We employed (1x1) period- 

icity, fixing the IML-GaN laterally and height to relax 
calculate the total energy. We were shown that the N- 
polarity structure was more stable than the Ga-polarity 
structure in Table The most stable for the IML- 
GaN on the 0-tcrminated ZnO(OOOl) interfacial struc- 
ture was N-polarity Wurtzite type interfacial structure 
shown in Figure ^b). The relative energy of the Ga- 
polarity Zincblend type interfacial structure shown in 
Fig. n^a) that was suggested from atomic diffusion was 
0.90eV. Therefore, it was shown that the GaN growth 
layer on the 0-terminated ZnO(OOOi) surface has the N- 
polarity Wurtzite interfacial structure. It was shown that 
both the Ga-polarity and N-polarity GaN possibilities 
of grown on the 0-terminated ZnO(OOOi) surface using 
the first principle calculation. We discuss the cause into 
which the polarity changes by paying attention to the mi- 
gration paths and the migration barrier energies. We are 
given the structure of the 0-terminated ZnO(OOOl) sur- 
face, we focus our attention on the relative stability of the 
four high-symmetry sites, namely, the original wurtzite 
and zincblende lattice site (L), the bridge site (B), the 
hollow and three neighbors site (H3), and the top of the 
second layer and four neighbors site (T4). 

For Ga adatoms the most stable adsorption site on the 
0-terminated ZnO(OOOl) surface is the H3 site. The en- 
ergetically lowest transition path is the B site and the T4 
site. The T4 site is significantly higher in energy of 0.49 
eV. Therefore, to hop from one T4 site to the next, the Ga 
adatom diffuses along the B site over the T4 site which 
is the transition site leading to the migration barrier en- 
ergy of 0.49 eV. The L site is unstable for adsorption of 
Ga adatom: the relative energy is 0.96 eV. The origin 
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FIG. 1: Atomic structure (top and side view) for the IML- 
GaN on the 0-terminated ZnO(OOOl) surface, (a) is the Ga- 
polarity zincblende interfacial structure, and (b) is the N- 
polarity wurtzite interfacial structure. 

of the energy difference comes from the fact that an iso- 
lated Ga atom has three bonds. At the H3 site, the Ga 
adatom uses all three bonds with three O atoms of the 
topmost layer. However, at the L site, the Ga adatom 
uses only one bond for adsorbation and the remaining 
two bonds are dangling-bonds. Thus, the crystal growth 
will be very difficult after constructing such a structure. 

For N adatoms the most stable adsorption site on the 
0-terminated ZnO(OOOl) surface is L site. The energet- 
ically lowest transition path is the B site. The B site is 
significantly higher in energy of 0.15 eV. Therefore, to 
hop from one L site to the next, the N adatom diffuses 
along the B site which is the transition site leading to 
the migration barrier energy of 0.15 eV. The origin of 
the energy difference comes from the fact that an iso- 
lated N atom has three bonds. At the H3 and T4 sites, 
the N adatom uses all three bonds with three O atoms 
of the topmost layer. However, the N-0 bond length is 
shorter than the Zn-0 bond length. 

Thus, the N adatom is bind to one O atom of the top- 
most layer. According to the calculation, an isolated N 
adatom is the most stable at the L site, and an isolated 
Ga adatom is unstable at the L site. Table HH show the 
migration barrier energies of the Ga and N adatoms were 
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TABLE II: The migration barrier energy values of the N on 
the Ga-terminated GaN(OOOl) surface and the O-terminated 
ZnO(OOOl) surface. The unit of the energy is eV. 



surface 


Migration barrier energy (eV) 




Ga adatom 


N adatom 


GaN(OOOl)" 
ZnO(OOOl)' 


0.4 


1.4 


0.49 


0.15 


"Rof. 5|. 
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compared with them on the Ga-terminated GaN(OOOl) 
smface that have been investigated theoreticaUy by Zywi- 
etz et al.5 It was shown that the migration barrier energy 
of the Ga adatom was similar to on the Ga-terminated 
GaN(OOOl) surface by 0.4eV. However, It was shown that 
the migration barrier energy of the N adatom was very 
smaher than on the Ga-terminated GaN(OOOl) surface 
by 1.4eV. Thus, the higher possibihty of the initial op- 
timized growth condition for the heteroepitaxy growth 
of GaN on the O-terminated ZnO(OOOl) surface is the 
Ga-polarity GaN crystal grown under the suitable low 
temperature condition. 

It was known to obtain the GaN crystal of the high 
quality when the GaN was grown under the Ga-rich 
growth condition. However, N adatom on the O- 
terminated ZnO(OOOl) surface diffuses compared with on 
the Ga-terminated GaN(OOOl) surface. Therefore, we 
have suggested that the adsorbed atomic specie of the 
first layer on the O-terminated ZnO(OOOi) surface is N- 
polarity because Ga adatom gives priority under the Ga- 
rich condition. Not Ga-rich growth condition but stoi- 
chiometric growth condition of the growth ratio is better. 
In general, when the GaN was grown on the Ga-polarity, 
it was known to obtain the GaN crystal of the high qual- 
ity. As for temperature of the initial growth condition, 
we were suggested that the stoichiometric growth condi- 
tion and the suitable low temperature condition be neces- 
sary compared with the typical growth temperature from 
those calculated results. 

We have experimentally investigated the polarity of 
GaN grown at low substrate temperatures under stoi- 
chiometric conditions by the use of PLD. Figure |21 shows 
a RHEED pattern for GaN grown at 300°C on a ZnO 
stepped and terraced substrate. GaN grown at this tem- 
perature shows a sharp streaky RHEED pattern, which 
indicates that it has high crystallinity and a flat sur- 
face. In addition, we have also observed RHEED oscilla- 
tion during the growth at 300°C, which implies that the 
growth proceeds in the layer-by-layer mode. It should be 
noted that RHEED image shown in Fig.[21has a x 2 recon- 
structed pattern, which is indicative of the Ga-polarity 
surface. Figures |31(a) and (b) show 1x1 fj,m^ AFM im- 
ages of the surfaces of GaN grown at RT before and 
after NaOH etching, respectively. A stepped and ter- 
raced structure can be seen in Fig.|3|(a), which is consis- 
tent with the layer-by-layer growth mode. One can see 



FIG. 2: A RHEED pattern for GaN grown at 300° C on a 
ZnO(QOOl) stepped and tellaced substrate. The x2 recon- 
structed surface is indicative of a Ga-polarity surface. The 
incidence of the electron beam is parallel to the GaN [1120] 
direction. 



that the surface of GaN remains atomically flat even af- 
ter the NaOH etching, which makes a concrete evidence 
for growth of the Ga-polarity GaN. The facts that GaN 
grown at substrate temperatures below 300° C under sto- 
ichiometric conditions has Ga-polarity and GaN grown 
at 600° C has mixed polarity as reported by Namkoong 
et al. 3 are quite consistent with our theoretical calcula- 
tions. 



FIG. 3: 1x1 fira^ AFM images of surfaces of GaN grown at 
RT (a) before and (b) after naOH etching. 



IV. SUMMARY 

In summary, we have investigated the stability of the 
IML-GaN on the O-polarity ZnO(OOOl) interface struc- 
ture using the first-principles calculation. We have found 
that the calculated results of the most stable of the 
IML-GaN on the O-polarity ZnO(OOOl) interface has 
the N-polarity. However, we have found that the calcu- 
lated results of the adatom diffusion on the O-terminated 
ZnO (0001) surface has the Ga-polarity. Therefore, the 
change in the polarity of GaN crystal grown on the 
O-terminated ZnO (0001) surface was important of the 
change of growth temperature. We have suggested that 
the optimized initial growth condition for the growth 
of the Ga-polarity GaN crystal on the O-terminated 
ZnO (0001) surface is under the stoichiometric growth 
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rate and the suitable low temperature condition. Experi- 
mental observations show that GaN grown on ZnO(OOOl) 
by PLD at substrate temperatures below 300° C has Ga- 
polarity, which is quite consistent with the theoretical 
calculations. 
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